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ABSTRACT 

Despite much theoretical and observational progress, there is no known firm upper limit to the 
masses of stars. Our understanding of the interplay between the immense radiation pressure pro- 
duced by massive stars in formation and the opacity of infalling material is subject to theoretical 
uncertainties, and many observational claims of "the most massive star" have failed the singularity 
test. LBV 1806—20 is a particularly luminous object, L~10® L©, for which some have claimed very 
high mass estimates (Minitiai>200 M©), based, in part, on its similarity to the Pistol Star. We present 
high-resolution near-infrared spectroscopy of LBV 1806—20, showing that it is possibly a binary sys- 
tem with components separated in velocity by ^70 kms~^. If correct, then this system is not the most 
massive star known, yet it is a massive binary system. We argue that a binary, or merged, system is 
more consistent with the ages of nearby stars in the LBV 1806—20 cluster. In addition, we find that 
the velocity of Vlsr=36 kms~^ is consistent with a distance of 11.8 kpc, a luminosity of 10^'^ L©, 
and a system mass of ~130 Mq. 

Subject headings: binaries: spectroscopic — stars: evolution — stars: early-type — stars: fundamental 
parameters — stars: individual (LBV1806-20) — infrared: stars 



1. INTRODUCTION 

Massive stars are key ingredients and probes of astro- 
physical phenomena on all size and distance scales, from 
individual star formation sites, such as Orion, to the 
early Universe during the age of reionization when the 
first stars were born. As ingredients, they control the 
dynamical and chemical evolution of their local environs 
and individual galaxies through their influence on the 
energetics and composition of the interstellar medium. 
They likely play an important role in the early evolution 
of the first galaxies, and there is evidence that they are 
the progenitors of the most energetic explosions in the 
Universe, seen as gamma ray bursts. As probes, they 
define the upper limits of the star formation process and 
their presence likely ends further formation of nearby 
lower mass stars. They are also prominent output prod- 
ucts of galactic mergers, starburst galaxies, and AGN. 

Despite their central role in these varied topics, there 
is no known firm upper limit to the maximum stellar 
mass. Such a basic quantity escapes both theory, because 
of the complex interplay between radiation pressure and 
opacity, and observation, because of incompleteness in 
surveying the Galaxy along the plane and the dearth of 
resolved starburst clusters with masses >10^ M©. The 
mos t promising cases ar e made recently bv iFieed (l2003f) 
and iWeidner fc Krouoal l)2004|) . who give observational 
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evidence for an upper mass cutoff of ^150 Mq in the 
Arches cluster and R136, respectively. Finding even a 
single star with a mass significantly beyond this limit 
would single-handedly change our understanding of the 
maximum mass a star may have. 

Of course, each "discovery of the most massive star" 
must face the crucial test of singularity. Some such 
claims have not passed this test, with subsequent obser- 
vations revealing that the obj ect is actually compos ed of 
multiple stars. For instance, iWeieelt et al.l l)1991l) find 
that R136 is actually a massive cluster of stars, rather 
than a supermassive single star with M^250-2,0 00 Mq 
(iFeitzinger. Schlosser. Schmidt-Kaler. fc Winkled 
[ 19801 IPanagJa. Tanzi. fc Tarenghil I1983|) . Likewise, 
IDam ineli et al.l lj2000|) claim that ri Car is binary, with 
present-day masses of ~70 Mq and ^^30 Mq, rather 
than a single star with M^lOO Mq. 

Some of the most massive sing l e star s known are near 
the Galactic Center. iFiger et all l)1998() identify the Pis- 
tol Star, in the Quintuplet Cluster, as one of the most 
massive known, with an initial mass of ~200 Mq, and 
show that the star is single based upon their Keck speckle 
data and spectra; the former reveal that the star is single 
down to a projected distance of 110 AU (14 mas), while 
the latter do not show an obviously composite spectrum. 
iFigcr ct al. ( 1999a) identify a near-twin to the Pistol 
Star, FMM362 less than 2' from the Pistol Star, and 
iGeballe et alJ lj2000(l demonstrate that it had the same 
spectroscopic and photometric characteristics as the Pis- 
tol Star, at the time of the observations. Presumably, 
F MM362 is nearl y as ma ssive as the Pistol Star. 

IKulkarni et alJ ljl995f) identify a luminous source 
near the soft gamma ray repeater, SGR 1806 —20. 
Ivan Kerkwiik. Kulkarni. Matthews, fc Neueebaueil 
(J1995J) determine that this source has characteristics 
similar to Luminous Blue Variables (LBVs), i.e. it has 
L^IG^ Lq and a K-band spectrum indicating a spectral 
type in the range 09-B2, leading them to propose that 
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it is an LBV-candidate. lEikenberrv et all l|2001l 12004(1 
estimate the properties of the star, finding a luminosity, 
and thus mass, that is at least compar able to that of 
the Pistol Star. lEikenberrv et alJ l|2004D obtained high 
resolution images at Palomar, showing that the star is 
single down to about 900 AU. Note that the putative 
binary comp onents of r] Car have a maximum separation 
of 30 AU l|Damineh et alJl2000() . So, it is possible 
that LBV 1806—20 is multiple, but it is certainly not 
a "cluster," in the traditional sense. Even the most 
compact cluster in the Galaxy, the Arches cluster, h as a 
half-hght radius of -40,000 AU (Fi ger et aUlToOObl) . 

In this Letter, we present high-resolution near-infrared 
spectra of LBV 1806—20 showing that the helium ab- 
sorption lines are double, suggesting that the star may 
be binary. 

2. OBSERVATIONS AND DATA REDUCTION 

The data were obtained at Ke ck on 22 June 2003, using 
NIRSPEC l|McLeanet aLllTool . LBV 1806-20 was ob- 
served at high resolution (R'^26,500=A/AAfwhm) from 
1.68 fiia to 2.25 fim in three grating settings, and at 
low resolution (R'^2,500) from 2.7 /im to 4.2 /j,m in two 
grating settings. We obtained two exposures per grat- 
ing setting, with the telescope nodded along the direc- 
tion of the slit inbetween exposures. The exposure times 
were 200 seconds, 100 seconds, and 30 seconds, for the 
high resolution short-wavelength, long-wavelength, and 
low resolution pointings, respectively. 

Quintuplet Star #3 (hereafter "Q3"), which is feature- 
less in this spectral region (Fiee r etaDllQQSL Figure 1) , 
was observed as a telluric standard ("Mon eti et alJll994f) . 
Arc lamps containing Ar, Ne, Kr, and Xe, were observed 
to set the wavelength scale. In addition, a continuum 
lamp was observed through an etalon filter in order to 
produce an accurate wavelength scale between arc lamp 
lines and sky lines (predominantly from OH). A quartz- 
tungsten-halogen lamp was observed to provide a "flat" 
image which was divided into the background-subtracted 
target images. 

We subtracted background flux, dark current, and 
residual bias by differencing the nod pair frames, and we 
divided the result by the flat fleld frame in order to com- 
pensate for non-uniform response. We extracted spectra 
of the objects and applied a wavelength solution based 
upon the locations of arc lines. Finally, we divided the 
target spectra by the spectrum for Q3 in order to remove 
telluric effects. Note that we omitted this last step in the 
case of the spectra near 1.7 /im because Q3 is too faint 
at such short wavelengths. The flnal extracted spectra 
are shown in Figures ^ 13 El a-nd 01 They have been 
individually shifted along the x-axis so that the spectral 
features are aligned with their vacuum wavelengths at 
rest. 

3. DISCUSSION 

3.1. Spectral Analysis 

The spectra of LBV 1806-20 and the Pistol Star are 
virtually indistinguishable, except for the morphology of 
the Fe II and He I lines. The latter are generated in the 
photospheres and winds of LBVs and arc invariably seen 
in absorption. In LBV 1806—20, the He I lines are obvi- 
ously double, with a velocity separation of 72±3 kms^^. 
The short-wavelength component of the 1.70 /im line has 
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Fig. 1.— Spectra of the Pistol Star and LBV 1806-20 near 
the 1.70 fim He I transition. The feature in the spectrum for 
LBV 1806—20 is double (marked by arrows), with the blueward 
component having a P Cygni signature. The broad feature on 
the left at 1.688 iim is attributable to an Fe II transition. The 
two features at 1.691 ^m and 1.698 fim are due to Si II. The 
weak absorption features near 1.6905 fim and 1.6960 fim are due 
to telluric absorption. 
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Fig. 2. — Spectra of massive stars near the 2.112/2.113 fim 
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The Pistol Star, IRS16NE. and LBV 1806-20 are 



tes IFiger e t al. 1998; Tamblvn et al. 199^ 
ik. K ulkarni. Matthews . & Neugebaucr 1995), and the 
are B-type supergiants I Figer et al. 19993 
iGenzel et aljil99M . All the spectra show a single He I doublet 
near 2.1125 fira, except that for LBV 1806—20, in which it is obvi- 
ously double (marked by arrows). The broad feature near 2.089 /^m 
is common to LBV spectra and corresponds to an Fe II transition. 
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Fig. 3.— Spectra of the Pistol Star and LBV 1806-20 near the 
Brackett— 7 line. The two absorption lines near 2.160 fira in the 
spectrum for LBV 1806—20 (marked by arrows) are due to two 
separate transitions of He I. 




LBV 1806-20 



0.0 t 




3.50 3.60 3.70 3.80 3.90 4.00 4.10 
Wavelength (urn) 

Fig. 4.— Spectra of the Pistol Star and LBV 1806-20 in the 3- 
4 fira region. Both spectra have a strong Brackett— a emission hne 
near 4.05 fiia, a Pfund— 7 (H I 6-5) emission line near 3.74 fim, and 
a set of emission lines in the Humphreys series (H I n-6). Hydro- 
genic transitions of He I can be seen in absorption just shortward 
of the Brackett— o lines. The remaining lines are mostly due to 
Fell. 
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Fig. 5. — Galactic rotation curve toward the line of sight 
of LBV 1806-20 from Brand & Blitz (1993). The cross-hatched 
region corresponds to our estimate of the systemic velocity of 
LBV 1806—20 in the LSR frame. The error band limits are 
±10 kms~^, where 6 kms~^ represents the measurement error, 
and 4 kms~ allows for mass ratios between 1:1 and 1:1.5. 



a P Cygni signature, indicating that the star responsible 
for the bulk of the emission lines is blue-shifted with re- 
spect to the other star. Indeed, none of the emission 
lines are double, indicating that the red-shifted compo- 
nent likely has less luminosity and/or a thinner wind. 
The velocity separation between the two minima in the 
components is 106 kms^^, larger than the separation of 
the components near 2.112 /ini, as would be expected 
given the P Cygni nature of the profile in the blue com- 
ponent of the 1.700 /im line. The He I absorption com- 
ponent near 2.160 /xm appears to be double (marked by 
arrows), with a separation of 99 kms^"'^; however, we 
believe the lines correspond to two expected He I lines 
with velocity seperation of 83 kms""'^ near that wave- 
length. Figure 0] does not appear to show any double 
absorption lines in the spectrum of LBV 1806—20, pre- 
sumably due to the relatively low spectral resolution of 
the data (70 kms""'^ corresponds roughly to 2 pixels in 
the low resolution mode). The absorption feature near 
4.04 fim is likely due to a blend composed of a He I 
doublet {5fFo-U^'D at 4.0377 ^lm and 5i^Fo-M^B at 
4.0409 /im). The He I absorption blend near 3.703 /im 
(5d3D-4p3Po) has a FWHM of 183 kms-\ compared to 
a FWHM of 152 kms~^ for the nearby Pfund-7 line at 
3.74 /im. 

The Fe II lines are broad in the spectra of both the 
Pistol star and LBV 1806—20, but they have multiple 
peaks in the latter. The 1.688 /xm and 2.089 /im lines 
have nearly identical morphology, each with three peaks 
in similar intensity ratios and relative positions in the 
spectrum of LBV 1806—20. The Gaussian shape of the 
Fe II lines in the spectrum for the Pistol star is likely 
caused by recombination closer to photosphere, where 
the wind velocity has not reached its constant terminal 
value. In the spectrum for LBV 1806—20, the iron likely 
recombines further out in the wind so that emission is 
dominated by gas at the terminal velocity, thus causing 
the flat top shape of the emission line. This explanation 
for the variation in the Fc II line shapes in the two stars 
suggests that LBV 1806—20 has slightly higher ionizing 
radiation, consistent with the apparent P-Cygni nature 
of the blue component of the He I lines in the spectrum 



for that star. 

Spectra of both stars have similar morphology in the 
H lines. However, the presence of deeper absorption fea- 
tures in the observed H lines of the LBV 1806—20, as well 
as their narrower nature, are consistent with underlying 
absorption H lines from the putative companion. 

We estimate the line-of-sight velocity for 
LBV 1806—20 by comparing the vacuum wavelength 
of the He I line near 2.1125787 /im with its measured 
positions at 2.1125065 /xm and 2.1130192 /im, as 
determined by fitting Gaussian profiles to the two 
absorption components. We adjust the velocity to 
the heliocentric frame by adding -1-0.2 kms~^ to the 
measured velocity to compensate for Earth motion at 
the time of observations. We add -1-10 ± 0.4 kms~^ 
to this value in order to a djust the velocity to th e 
local standard of rest frame l)Dehnen fc BinnevI 11998(1 . 
Assuming a roughly equal mass binary, we average the 
velocities of the two absorption lines, to find a systemic 
velocity of Vlsr=36.4 kms"-'^; note that a change in this 
assumption, to a mass ratio of 1:1.5, would only induce 
a change in our estimate of ±4 kms~^. Our estimate is 
similar to those we obtain for three nearby hot stars, for 
which Vavc,LSR=35.0±5.7 kms^^ (see Figure HI). These 
veloci ties contrast with the value in lEikenberrv et all 
l|200^ of 10±20kms-i. 

The bulk velocity of the LBV 1806-20 cluster can 
be inter preted in the c o ntext of the Galactic rotation 
curve in iBrand fc Blitd 1)1993(1 in order to estimate its 
distance. We find a distance of 4.1 kpc or 11.8 kpc 
(see Figure E)). Only the high solution is consistent 
with a location beyond an NH3 cloud at a distance of 
5.7 kpc (Brand fc Blitz 1993) . This estimate is c onsid- 
erably less than the estimate in lEikenberrv et aD l(2004fl 
of 15.lli;^ kpc. 

3.2. Mass of LBV 1806-20 

For the far distance estima te (d ^ ll.8 kpc), and 
the photometry in Eikcnbcrrv et al.l 1(20041) . we find 
Mx=— 9.5±0.1. It is difhcult to convert this value into 
a luminosity estimate because the range of tempera- 
tures, and, therefore the K-band bolometric correction. 
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varies widely for LB Vs l|Blum. Depov. fc SellgreiJfl995D . 
However, if we us e the BCk for the Pistol Star from 
iFiger et al] l)1998|) . then we arrive at a luminosity of 
10^'^ Lp,. Using Langer's stellar evolution models in 
IFiger et alJ l)1998D . we find an initial mass '^130 Mq in 
the case that the object is a single star. 

3.3. LBV 1806-20 as a Massive Binary 

We suggest that the object is a massive binary, based 
upon the double lines in its spectrum. For a 1:1 mass 
system, the components have initial masses of 65 Mq, 
assuming a mass-luminosity relation of McxL", where 
alpha'^1. Our present analysis is fairly insensitive to the 
mass ratio, a topic that we will revisit in a future study. 
By comparison, the initial m asses of the components in 
r? Car are 120 Mfn and 40 Mrr-. (iDamineh. Conti. k LopesI 
Il997t iPittard fc CorcoranI 120021) . Another massive bi- 
nary system, WR20a, contains roughly equal mass com- 
ponents, each haviii K M--70-80 Mq ( Rauw et al. 2004; 
IBona,noset 711 12004;) . 

LBV 1806—20 is similar, in appearance, to a growing 
group of objects, i.e. the Pistol Star, FMM362, 77 Car, 
and others. Whether these objects are single or multiple, 
they all presumably have at least one evolved component 
that is massive (^50 Mq) and some are surrounded by 
circumstellar material produced by massive eruptions. 

3.4. The Formation of the Most Massive Stars 

Theoretical arguments suggest that stars with ini- 
tial masses ^150 Mq have hydrogen-burning (main se- 
quence) lifetimes of about 2 Myr (|Bond. Arnett. fc Card 
11984') ■ and total lifetimes no more than ~3 Myr 
(jMevnet et 11111994^ . If LBV 1806-20 is single, and it 
had an initial mass in this range, then it would necessar- 
ily be much younger than nearby Wolf-Ra yet stars in the 
same cluster that a rc at least 3 Myr old IjMevnet et al.l 
I1994J) . lEikenberrv et al. (2004'! suggest that a supernova 
from a first generation of stars could have triggered the 
formation of LBV 1806—20, and perhaps other nearby 
stars. While this mechanism is often invoked to explain 
mutli-generation star formation, it is not the only expla- 
nation in this case. Note that the LBVs in the Quin- 



tuplet cluster should not live longer than ^^3 Myr old, 
whereas the other cluster members are 4 Myr old; how- 
ever, in this case, the original natal cloud must have been 
strongly sheared by the tidal field in the Galactic center 
on time scales much less than 2 Myr. A more general 
explanation for the age discrepancy in both the Quintu- 
plet and in the LBV 1806—20 cluster might be that the 
most massive stars are binary systems of less massive 
stars (which can live longer), or that they were binary 
systems that recently merged (jFiger fc Kimll200^ . 

This proposition also resolves the problem of explain- 
ing why 200 Mq stars are observed in the Quintuplet 
cluster, whereas the upper mass cutoff to the initial mass 
funct ion in the nearby Arches cluster is ~150 M© (Figen 
|2003)- In fact, if single stars with masses of ^100 Mq 
exist, then where are the binary systems with such stars? 
Certainly, some of the most massive stars exist in binary 
systems, as evidenced by the substantial binary fraction 
( 25%) in the WR catalogue ign dcr Hucht 2001). Per- 
haps we are witnessing the ancestors of such systems in 
these LBVs. 

The evidence in this paper suggests that LBV 1806—20 
may be binary, but we will continue to explore the 
possibility that the double-lined spectra could have 
an alternate explanation. For instance, perhaps the 
profile is a composite of an emission feature, gen- 
erated by a disk or wind, at the systemic veloc- 
ity and superposed on a photospheric absorption fea- 
ture. In addition, we will refine the temperature, lu- 
minosity, and mass ratio estimates through a more de- 
tai lcd analysis, coupled with wind / atmosphere mod- 
els [Najarro^^^iger^^^^dntzkj l|2004fl . In addition, we 
have begun a spectroscopic monitoring program to de- 
tect changes that might indicate a binary nature for 
LBV 1806-20. 
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